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ABSTRACT 

We report follow-up XMM-Newton and ground-based optical observations of the unusual X-ray 
binary SDSS J102347.67+003841.2 (=FIRST J102347.6+003841), and a new candidate intermediate 
polar found in the Sloan Digital Sky Survey: SDSS J093249.57+472523.0. SDSS J1023 was observed in 
its low-state, with similar magnitude/color (V = 17.4 and B = 17.9), and smooth orbital modulation 
as seen in most previous observations. We further refine the ephemeris (for photometric minimum) to: 
HJD(TT) min = 2453081.8546(3) + E * 0.198094(l)d. It is easily detected in X-rays at an unabsorbed 
flux (0.01-10.0 keV) of 5 x 10~ 13 ergcm~ 2 s _1 . Fitting a variety of models we find that: (i) either a hot 
(kT > 15 keV) optically thin plasma emission model (bremsstrahlung or MEKAL) or a simple power 
law can provide adequate fits to the data; (ii) these models prefer a low column density ~10 19 cm~ 2 ; 
(iii) a neutron star atmosphere plus power law model (as found for quiescent low-mass X-ray binaries) 
can also produce a good fit (for plausible distances), though only for a much higher column « 4 x 10 20 
cm~ 2 and a very cool atmosphere (kT < 50 eV). These results support the case that SDSS J1023 is a 
transient LMXB, and indeed places it in the subclass of such systems whose quiescent X-ray emission 
is dominated by a hard power law component. Our optical photometry of SDSS J0932 reveals that it 
is an high inclination eclipsing system. From our two epochs of data, and 7 eclipse times, we are able 
to derive a best fit ephemeris for minimum light: HJD(TT) min = 2453122.2324(1) + £*0.0661618(4)d, 
although aliases, with one cycle count different between epochs, are acceptable. The X-ray spectrum is 
well fit by either a hard bremsstrahlung or power law, with a partial covering absorption model, with 
a high covering fraction ^0.9 and column w 10 2 



i-!3 cm 2 _ c om bined with its optical characteristics 



high excitation emission lines, and brightness, yielding a large Fx/F opt ratio — this highly absorbed 
X-ray spectrum argues that SDSS J0932 is a strong IP candidate. However, only more extensive optical 
photometry and a detection of its spin or spin-orbit beat frequency can confirm this classification. 

Subject headings: individual: (SDSS J093249. 57+472523.0, SDSS J102347.67+003841.2) — novae, 
cataclysmic variables — stars: magnetic — X-rays: stars 



1. INTRODUCTION 

Cataclysmic variables (CVs) and low-mass X-ray bi- 
naries (LMXBs) are two broad classes of interacting 
binary star, whose emission, from X-ray through opti- 
cal/IR (and even radio), is powered by accretion onto a 
compact object (white dwarf or neutron star/black hole 
respectively) . There are good reviews of the various sub- 

i Some of the results presented here were obtained with the 
MMT Observatory, a facility operated jointly by The University of 
Arizona and the Smithsonian Institution. 
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classes and general characteristics in lWarneil l)1995(K and 
iLewin et al.l l)1995|) respectively. When actively accret- 
ing, both classes exhibit (broad) emission line spectra, 
with very blue continua (from the accretion disk), but 
are normally distinguished by their relative X-ray out- 
put; their X-ray luminosity roughly scales as M*/i? 2 , i.e. 
a factor of 10 4 — 10 5 . However, in many cases the accre- 
tion rate onto the compact object is highly variable, with 
both CVs and LMXBs exhibiting a variety of brightness 
states. 

For the past three years, we have used XMM-Newton 
to study and confirm potential mag netic CVs identi- 
fied in the Sloan Digital Sky Survey ijYork et alJ 120001 
SDSS). Two SDSS sources stood out based on the un- 
usual strength of their He emission lines, as well as on 
the double-peaked character of all their lines. These 
are SDSS .1102347.67+00384 1.2 (ISzkodv et a.l.ll2003T) and 
SDSS J093249. 57+472523.0 llSzkodv et al .1120041) . 

SDSS J102347.67+003841.2 = FIRST J102347.6 
+003841 (hereafter SDSS J1023) was identified as a CV, 
from its radio em ission in the FIRST radio survey 
ijBond et alJ 120021 , the first CV to be so discovered) , 
and independently from its color and emission line spec- 
trum in the SDSS. In both cases, a highly magnetic white 
dwarf was suggested to account for the high-excitation 
optical emission lines and detectable radio emission. The 
nature of the lines appeared more consistent with an in- 
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termediate polar (IP), since they lacked both the large 
equivalent widths and narrow components usually seen 
in polars. Significant changes in state have been pos- 
tulated to account for: (i) very different opt ical photo- 
metric behav i our s een by iBond et alJ ((2002) and later 
iWoudt et alJ l)2004fl . a change from pure flickering be- 
haviour to a smooth periodic modulation on a period of 
4.75 hr; (ii) significant color changes between the SDSS 
photometry (when it appeared red) to the later SDSS 
spectroscopy, sh owing a very blue continuum. 

Most recently, T horstensen &: A rmstrong ( 2005) (here- 
after TA05) presented the results of a photometric and 
spectroscopic campaign, which have called into ques- 
tion whether SDSS J1023 is a CV, or a LMXB instead. 
The low-state spectrum reveals a mid-G star, unlike 
either previous spectral observation. Multi-color pho- 
tometry showed a s mooth orbital modulation similar to 
IWoudt et alJ l)2004[) . with color changes consistent with 
the changing aspect of the heated face of the secondary, 
but even at maximum light no emission lines are present 
in the spectrum. Modelling of the light curves, together 
with the radial velocity constraints, indicates that a high 
primary mass (greater than the Chandrasekhar mass) is 
needed, thus making SDSS J1023 a LMXB. We obtained 
X-ray, circular polarization and other optical data to fur- 
ther elucidate the nature of this source. 

For SDSS J0932 we also obtained contemporaneous op- 
tical light curves, spectra and X-ray/optical data from 
XMM-Newton in an attempt to discern its nature and 
possible relation to SDSS J1023. 

2. OBSERVATIONS 

A summary of the X-ray and optical observations is 
presented in Tabled XMM-Newton possesses three X- 
ray telescopes, bac ked by the two MOS and one pn 
ijTurner et al.ll2001[) CCD came ras, and an optical mon- 
itor fOM. IStriider et al.ll2001|k hence parallel data are 
obtained by all detectors. The two Re flection Grating 
Spectrographs ijden Herder et afl 1200 lj) are arrayed in 
the optical path of the MOS detectors; however, for nei- 
ther target were there sufficient counts to provide useful 
signal in the resulting dispersed spectra. In the direct 
(spectro)imaging EPIC detectors count rates were suf- 
ficient for extraction of low-resolution spectra and light 
curves. We followed the standard protocol as given by 
the Vilspa XMM-Newton Science Analysis System (SAS) 
site 9 and the ABC guide 10 from the US GOF. In both 
cases we used calibration files current to 2005 March 16 
and the SAS v6.1.0, and as precaution we used SAS-tools 
to produce new event list files from the Observation Data 
Files incorporating the latest calibration updates. 

To check on variations in the X-ray background we cre- 
ated light curves for each entire detector in the 10-15 keV 
range. No hard X-ray flaring events were present in the 
data, and for SDSS J1023 we simply used the standard 
good time intervals (GTIs). The case of SDSS J0932 
was more complex. Due to its low count rate, the effects 
of background are relatively important. Examination of 
the background spectra showed the usual steep rise at low 
energies (due to detector noise) and a number of promi- 



9 Available from http://xmm.vilspa.esa.es/external/ 
xmm_sw_cal/sas.shtml 

10 http:/ /heasarc. gsfc.nasa.gov/docs/xmm/abc/abc. html 



nent metal lines, which are due to fluorescence of satel- 
lite materials bombarded by soft protons, especially in 
the pn. For the pn it is possible to make some improve- 
ments to reduce this detector background; we followed 
the guidelines in SAS handbook section 4.3.2.1 "Improv- 
ing the quality of EPIC pn data: epreject." Next, we cre- 
ated background light curves to examine any variations 
at these specific energies, and identified one interval of 
especially high flux, and created new GTIs to exclude it. 

In addition to the time selection, we screened the event 
lists using the standard canned expressions and restricted 
energies to the 0.1-10 keV range. Owing to the very dif- 
ferent count rates, we used different extraction radii for 
each source. For SDSS J0932 we used circular source 
apertures of radius 320 pixel (MOS) and 360 pixel (pn), 
in each case enclosing ~ 70% of the energy, chosen to 
maximize S/N. For SDSS J1023 we increased the radii 
to 480 and 640 pix respectively, now enclosing ~ 80% of 
the energy. For the background, we selected annuli on 
the same chip for the MOS, and adjacent rectangular re- 
gions at similar detector Y locations to the target for the 
pn, excluding any point sources. As advised we selected 
single to quadruple events for the MOS (pattern < 12), 
but only singles from the pn (pattern = 0). 

Using the same extraction regions for source and back- 
ground we extracted their light curves with SAS task 
evselect. In this instance (to maximise counts) we in- 
cluded both singles and doubles from the pn (pattern< 
4) . Finally, background subtraction (scaled by extraction 
area) and correction of the time stamps to the heliocenter 
(times in HJD(TT) 11 ) were undertaken with FTOOLS 12 . 
We obtained simultaneous B-band light curves from the 
OM. Again, we performed our own extractions using 
omf chain starting from the ODF event list files, setting 
our binning times at this earliest stage. We converted 
the count rates into magnitudes according to the results 
from interactive photometry with the omsource task. 

SDSS J0932 was also observed by the US Naval Obser- 
vatory Flagstaff Station (NOFS) lm telescope on 2005 
April 11 and May 12. The latter observation covers the 
3 hr window immediately preceding the XMM-Newton 
observation (and OM optical coverage). To maximise 
signal-to-noise these data were taken in white light, but 
an approximate zeropoint for this differential photom- 
etry to Johnson ~ U-band was made possible through 
calibration of the field from all-sky photometry including 
Landolt standards observed at NOFS. Hence, the light 
curves are labelled V magnitudes. A night after its X- 
ray observations (2005 May 13) a similar time series was 
taken on SDSS J1023. Circular spectropolarimetry and 
spectrophotometry was also obtained for SDSS J1023 
with the CCD Spectropolarimeter SPOL (Schmidt et al. 
1992) in 2002 May on the 2.3 m Kitt Peak Bok telescope 
and on the 6.5 m MMT in 2004 February, as indicated 
in Table 1. At both epochs, the object displayed the 
spectrum of a G star, with apparent magnitudes of 17.9 
and 17.5, respectively. The 20 min observation sequences 
revealed no significantly polarized features, and respec- 

11 The tools actually yield barycentric Julian Date in the 
barycentric dynamical time system, BJD(TB). However, the offset 
to heliocentric Julian Date in the geocentric (terrestrial) dynamical 
time system (HJD(TT)) is less than ~3 s at any given time, fine 
for our purposes here. 

12 http:/ /heasarc. gsfc.nasa.gov/lhcasoft/ftools/ 
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TABLE 1 

Observation Summary 





tit riatp 


Obs 


UT Time 




( Vtrvi in iQTl^~G 


1023 


2002 May 11 


Bok:SPOL 


03:54 - 


04:17 


circ. pol.< 0.5% 


spectropolarimetry, 1200s 




2004 Feb 16 


MMT:SPOL 


08:40 - 


09:02 


circ. pol.< 0.03% 


spectropolarimetry, 1200s 




2004 May 12 


XMM- Newton: 














EPIC-pn 


10:44 - 


14:37 


0.08 cts s" 1 


12570s live time b 






EPIC-MOS1/2 


10:21 - 


14:42 


0.03 cts s" 1 


15430s live time 






OM 


10:30 - 


14:14 


B = 17.9 


12899s duration 




2004 May 13 


NOFS 


03:12 - 


06:55 


~ V = 17.3 - 17.6 


open filter photometry 




2004 May 23 


APO: DIS 


05:40 - 


06:05 




spectrum (poor fluxing) 


0932 


2004 Apr 11 


NOFS 


03:53 - 


08:48 


~ V = 18.5 - 19.0 


open filter photometry 




2004 May 12 


NOFS 


03:08 - 


06:41 


~ V = 18.5 - 19.2 


open filter photometry 




2004 May 12 


XMM-Newton: 














EPIC-pn 


07:08 - 


08:57 


0.008 cts s' 1 


5698s live time 






EPIC-MOS1/2 


06:46 - 


02:00 


0.003 cts s _1 


7953/7984s live time 






OM 


06:54 - 


09:04 


B = 19.1 


7478s duration 




2004 May 23 


APO: DIS 


05:08 - 


05:33 


~ V = 19.3 


spectrum 



a The open filter photometry from NOFS has an estimated equivalent V zero-point, while for spectra 
the flux density at ~5500A is used. The XMM-Newton count rates are average values for each 
observation for a single detector. 

b The live time of the X-ray CCD detectors refers to the sum of the good-time intervals, less any 
dead time. It is typically much less than the difference of observation start and stop times. 



tive upper limits are v = V/I = 0.5% and 0.03% when 
summed over the spectrum. 

We also observed both SDSS J0932 and SDSS J1023 
on 2005 May 23 with the double-imaging spectrograph 
(DIS) with a resolution of about 2A on the Apache Point 
Observatory (APO) 3.5m (see fig.Q). These spectra were 
reduced and flux-calibrated using standard IRAF 13 rou- 
tines. 

3. SDSS J1023 

3.1. Light Curves 

Figure [21 shows that when we combine our two opti- 
cal light curves (phase-folded according to our refined 
ephemeris), we have full phase coverage. Again we see 
the same kind of smoot h (but non-sinusoid al) modula- 
tion as observed by both lWoudt et alJ l)2004|) and TA05. 
Both our optical curves indicate a broad (~0.3 in phase) 
constant flux interval at maximum light, followed by a 
smooth rounded minimum. Excellent agreement is seen 
for the phasing of minimum light (posited as superior 
conjunction of the donor); 574 cycles have elapsed be- 
tween the TA05 ephemeris and our epoch of observation, 
giving a phase uncertainty of merely 0.006. Fitting a 
Gaussian to the minimum in the 2004 May 13 NOFS 
lightcurve we determine Tq = 2453138.7077(3), giving 
P' = A(T )/574 = 113.7062(6)/574 = 0.198094(l)d, and 
a refined ephemeris of: 

HJD(TT) min = 2453081.8546(3) + E * 0.198094(l)d 

Note our ephemeris is given in Terrestrial Time (= UT 
+ 64.184s at the current epoch). 

13 IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy 
(AURA) Inc., under cooperative agreement with the National Sci- 
ence Foundation 

13 We use Terrestrial Time throughout this paper, since it the 
system best suited for ephemerides, unlike UT which is adjusted 
by leap seconds as needed to match the earth's rotation 



We constructed the X-ray light curve (Fig. [5} by sum- 
ming the background-subtracted light curves from all 
three EPIC cameras, with 300s bins. The X-ray light 
curve does not show a similar smooth modulation, but 
instead is quite irregular. Indeed, the amplitude of its 
variability (peak-peak) amounts to about 60%, greater 
than the 40% in B or 30% in white light. In view of the 
IP nature hypothesis, we note that it is possible to fit 
a sinusoid plus first harmonic model to our data, with a 
period of 0.069±0.003d (98 min), or0.3xP orb . As to the 
significance of this signal, th e peak in the Lomb-Scargle 
periodogram l)Scarglej lT982) is at the 96.5% confidence 
level, but we caution that with merely two cycles of cov- 
erage its reality remains debatable. Moreover, only three 
IPs have similarly large spin-orbit period ratios, V1025 
Cen, DW Cnc and EX Hya, and all these are IPs be- 
low the period gap, with periods ~1.5 hrs. If this pu- 
tative spin period were to be confirmed, it would place 
SDSS J 1023 in a unique position in the P or b — P S pin dia- 
gram. 

3.2. X-ray Spectral Fitting 

Following the most recent guidelines on the low-energy 
EPIC calibrations, we restricted our fitting to > 0.2keV 
for the MOS and > 0.15kcV for the pn, whilst still tak- 
ing care to check the effects of different choices of these 
limits. After initial inspection of the data, showing a 
high energy tail, we used all good data below 15 keV for 
this relatively bright source. There were sufficient counts 
to bin the data at >20 counts/bin, and use x 2 statis- 
tics to find the best fits to the background subtracted 
source spectra. We also double-checked our results us- 
ing >5 counts/bin, and fitting simultaneous source and 
background spectra using Cash statistics; in all cases the 
results were consistent with the x 2 approach. 

A variety of models were fit to the data, starting with 
simple single component plus interstellar absorbing col- 
umn (see Table. 0). In general the emission of CVs can be 
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Fig. 1.— APO spectra of SDSS J0932 (left panels) and SDSS J1023 (right panels). SDSS J0932 was observed in a similar state to 
that seen during its SDSS spectroscopy; it exhibits strong He lines indicative of high-ionization states seen in magnetic CVs. In contrast, 
SDSS J1023 was observed in its low state, showing neither blue light from an accretion disc nor any emission lines. The flux scale is in 
units of flux density 10 — 16 erg cm -2 s" 1 A 



approximated with a thermal plasma. We performed fits 
setting Nh — 4.22 x 10 20 cm~ 2 , the maximum Galactic 
column as derived from dust maps (using HEAsoft tool 
nH) or a fiducial Nh = 1 X 10 19 cm~ 2 , appropriate for 
such a high Galactic latitude source, and lastly with this 
parameter free. A high temperature (> 30keV) thermal 
plasma (either bremsstrahlung, or bremsstrahlung plus 
lines- mekal) can produce an acceptable fit (see fig. |3Jl 
but only with the lower column option (leaving the col- 
umn free yielded only an upper limit of Njj = 4.4 x 10 19 
cm , but was consistent with zero). The best fit was 
a simple power law with \t ~ 1-0 for either fixed col- 
umn, and photon index of 1.3—1.5 (also see fig. EJ- In 
this case, the data could roughly constrain the column, 
when freed, to Nh = 9 ± 7 x 10 19 cm~ 2 . A power law 
fit is suggestive of that found f or the hard componen t of 
LMXBs in quiescence (see e.g. lCampana et aI1ll998j) . 

Moving onto more complex models, we exam- 
ined two temperature and multi-temperature mekal 
(XSPEC:CEMEKl) models-used to fit IP spectra, and 
lastly a soft blackbody + mekal as for a polar. A two 
temperature fit was possible, but in the case of the low 
column the higher temperature mekal pegged at its up- 
per limit of 80 keV, whilst the cool plasma settled at 
100-300 eV. The data were not suitable to constrain the 
CEMEKL model at all. Adding in a soft blackbody to a 
MEKAL fit has a similar effect as the two temperature 
model. We found well-constrained kT-QB = 130 and 220 



eV respectively for the high and low column cases, but in 
both cases ^Imekal = 80keV (pegged). This compares 
to a typical polar fit with ^30eV and 30 keV for these 
components. 

On the other hand if SDSS J1023 is a LMXB, then we 
might expect a good fit for the typical two component 
neutron star atmosphere (NSA) plus power law model. 
Fixing both the neutron star radius and mass at their 
canonical 10 km and 1.4 M© values, the only free param- 
eters left in the NSA model are the distance and tem- 
perature. TA05 provide a constraint on the distance to 
SDSS J1023 depending on the mass of the secondary. We 
investigated fits for a range of masses/distances from 0.1 
M (1.O kpc) to 1.2 M Q (2. 3 kpc). Statistically acceptable 
fits are found for all distances, although in all cases only 
for a higher column > 10 20 cm -2 . However, especially 
for d — 1 kpc, but also for the longer distances, the ef- 
fective temperature of the NSA is low, and in all cases 
it contributes merely ~10% of the flux, even including 
energies down to 0.01 keV. If one performs an F-test to 
compare the best fit PL (which has a low column) to the 
best fit PL+NSA (which requires a higher column), the 
low value indicates that there is no need for the inclusion 
of the soft component to model the data; but of course 
this all depends on the actual column. 

In general, the quality of the fits is fairly similar and 
each yields an unabsorbed 0.01-10 keV X-ray luminos- 
ity of 5 x 10 -13 erg cm -2 s _1 . Although two-component 
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TABLE 2 

X-ray Spectral Fits for SDSS J1023 



Model 




reduced 

X 2 x 


10 20 cm- 2 


fcT a 
(eV) 


kT or T 
(keV) 


Flux b 




bremss 






1 


45 


4 


22 c 




15±* 


4.9 ±0 


4 








1 


06 





1 




63+f 


5.2 ±0 


8 


MEKAL 






1 


38 


4 


22 




12+1 


5.1±0 


5 








1 


02 





1 




Jy -12 


5.1±0 


6 


PL 






1 


16 


4 


22 




1.47 ±0.04 


5.5 ±0 


3 








1 


02 





1 




1.27 ± 0.03 


5.3 ± 


3 








1 


02 





q -t- n 7 




l.Ol It u.uo 


^ ^ -t- n 
0.0 IE u 





2 T MEKAL 






1 


09 


1 


22 


98+4 


23± 7 4 


11±4 













99 





1 


300l 2 * 


80 +pegged 


c: 9+O.2 
°' z -0.4 




BB + MEKAL 






1 


00 


4 


22 


130 ± 20 


80 +pcggod 


5.5 ±0 


3 











98 





1 


220 ± 40 


gg+pcggod 


5.0 ±0 


3 


PL + NSA: 


1 


kpc 


1 


02 


4 


22 


33 ± 1 


1.32 ±0.05 


6.0 ±0 


6 






1 


02 





o+3.2 
8 -0.7 


< 21 


i Q 9 +0.04 
1 "- >z -0.05 


5.2 ±0 


5 




2 


kpc 


1 


00 


1 


22 


45 ±2 


i 97 +0.05 
1 - z '-0.06 


5.8 ±0 


6 








1 


01 


4 


q+1-8 
°-2.1 


45 ±5 


1 24+ ' 07 
1,z -0.04 


5.9 ±0 


6 




2.3 


kpc 





99 


4 


22 


48 ±2 


1 26+ 05 
1 ' ZD -0.06 


5.8 ±0 


6 








1 


02 


4 


4 +l.6 
*-1.4 


48 ±6 


i o 9 +0. 06 
-0.04 


5.8 ±0 


6 



a For neutron star atmosphere (NSA) model the effective temperature (for a 
distant observer) is quoted, i.e. kT c s = g r x T c g where g r = (1 — 2.952 x 
Af na /_R ns ) 0,5 is the gravitational redshift parameter. 

b Unabsorbed flux in the 0.01-10 keV range, in units of 10 -13 erg cm -2 s _1 . 

c Model parameters without quoted uncertainties have been fixed at the value 
given. 



models yield slightly better fits, in most cases consid- 
eration of the parameter values indicate problems: ex- 
treme plasma temperatures or unusually low/high tem- 
peratures for the neutron star (either NSA or blackbody) . 
We defer further discussion to § [5] 

4. SDSS J0932 

4.1. Light Curves 

Owing to the low X-ray count rate of this source it was 
not possible to construct light curves with any useful time 
resolution. However, at B ~ 19 the OM light curve could 
usefully be binned to 160s, similar to the time-resolution 
of the NOFS data. We present the three light curves from 
2004 April and May in figure 0] Consistent with the 
double-peaked optical spectra, this source is clearly at a 
high-inclination to our line of sight, exhibiting deep, 2 
magnitude eclipses. We have undertaken anO-C anal- 
ysis using the 7 eclipse times spanning roughly a month. 
From the 4 May points alone we can determine the pe- 
riod to within 0.2%. Unfortunately, this is not quite 
high enough precision to facilitate an unambiguous cy- 
cle count back to the April epoch. The best fit solution 
yields: 

HJD(TT) min = 2453122.2324(1) + E * 0.0661618(4)d. 

In table 01 and fig. El we present our full 
O — C results for the 3 candidate solutions. Clearly, addi- 
tional photometry will be needed to confirm the correct 
one. 

SDSS J0932 exhibits variability in addition to its 
eclipses, which is notably different between the two 
epochs. In 2005 April there appears a long term trend, a 



slow rise and decline, with no obviously repeatable mod- 
ulation on the orbital period. In contrast, in 2005 May, 
there is no long timescale change in flux, but there is a 
clear periodic smooth variation. Close inspection of the 
light curves in the lower panel of fig. 0] reveals a possi- 
ble trend in the appearance of the pre-eclipse hump: it 
appears slightly later in phase each cycle. We fit a si- 
nusoid plus first harmonic model to all the non-eclipse 
datapoints, this gives a best fit period of 0.0688±°;™^ 
d, which is 3.8% longer than the orbital period, hence 
nominally discrepant at the 4cr level. Again with merely 
four cycles of data, this period is far from conclusive. In 
fig. El we show the NOFS data folded on both the orbital 
period and this longer one. The obvious explanation for 
such a period a few perc ent longer than the orbital is 
the super hump effect (see lPattersonlll99"8l for a detailed 
discussion of the supcrhump phenomenon). 

4.2. X-ray Spectra 

Even after our application of epreject, the back- 
ground is estimated to contribute 45% of the pn flux 
within our source aperture (only < 15% for the MOS 
detectors). With 50 and 55 counts total in the pn and 
MOS detectors respectively, the use of smaller 5cts/bin 
binning and Cash statistics is the best approach. First, 
we examined and fit the background spectra. For this 
we could bin to 20cts/bin and still have good spectral 
resolution, and we simply aimed to find the best empir- 
ical model to represent these contributions. Our final 
model consists of a twice-broken power law (i.e. 3 dif- 
ferent indices), plus Gaussian lines to account for the 
various metal emission lines. We performed a joint fit to 
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TABLE 3 

Optical eclipse timings and O — C results for SDSS J0932. 



HJD(TT) 


o-(HJD) 


P = 


0.0663035(4) d 


P = 


0.0661618(4)d 


P = 


0.0660206(5)d 


2453000 + 


(xlO- 4 ) 


Cycle 


O-C 


Cycle 


O-C 


Cycle 


O-C 






Count 


(cyclesx 10~ 5 ) 


Count 


(cyclesx 10 -5 ) 


Count 


(cyclesx 10 -5 ) 


106.68445 


4.5 


1 


3.3 ±6.8 


1 


0.0 ±6.8 


1 


-3.2 ±6.8 


106.75044 


2.0 


2 


-1.4 ±3.0 


2 


-2.5 ±3.0 


2 


-3.7 ± 3.0 


106.81690 


2.8 


3 


1.0 ±4.0 


3 


2.0 ±4.2 


3 


3.0 ±4.2 


137.64816 


1.6 


468 


2.9 ±2.4 


469 


0.5 ± 2.4 


170 


-1.9 ±2.4 


137.71423 


1.1 


469 


-0.5 ±2.1 


470 


0.8 ±2.1 


171 


-1.0 ±2.1 


137.78013 


2.3 


470 


-6.6 ±3.5 


471 


-4.7 ±3.5 


472 


-2.8 ±3.5 


137.8470 


3 


471 


1.4 ±4.5 


472 


5.5 ±4.5 


473 


9.6 ±4.5 



TABLE 4 

X-ray Spectral Fits for SDSS J0932 



Model 


Goodness a 


N H 


Partial 


covering 


kT (keV) 


Flux b 




of Fit 


(cm -2 ) 


N H (cm -2 ) 


frac. 


or r 




bremss 


70% 


&t\ X 10 22 






200 (pegged) 


2.4 ±0.8 




56% 


1.39 c x 10 20 


10^3 X 10 22 


0.95 ± 0.02 


200 (pegged) 


1.5 ±0.7 




61% 


1.39 x 10 20 


lOt^ x 1q22 


0.95 ± 0.02 


30 


1 5+ 16 


PL 


76% 


3 ± 3 x 10 22 






0.5 ±0.7 


d 




59% 


1.39 x 10 20 


9 ± 5 x 10 22 


0.9 ±0.1 


0.8 ±0.6 


1.7±0.6 



a As the fitting utilized Cash statistics, a Monte Carlo method was used to find the 
percentage of simulated spectra based on the parameter space of the model fit had 
C-statistic values less than than of the fit to the data. A good fit should have a value 
around 50%. 

b Unabsorbed flux in the 0.01-10 keV range, in units of 10~ 13 erg cm -2 s _1 
c Model parameters without quoted uncertainties have been fixed at the value given. 
d The singly absorbed power law is such an unphysical fit that flux estimates were 
meaningless. 



MOS background spectra, and independently for the pn. 
Next, we used this model (with relative normalisations 
scaled according to extraction area) as the fixed back- 
ground component in our source±background fitting. 

With our poor statistics we concentrated on find- 
ing the simplest model fit to the data, starting with a 
bremsstrahlung plus a single absorption (i.e. for Galac- 
tic) column. First, we found that there is significant 
column in excess of that estimated for the line-of-sight, 
our best fit requiring Nu = 5.6l 2 '3 x fO 22 cm -2 , as com- 
pared to the dust-map estimate (from the nH tool) of 
Nh = 1.39 x fO 20 cm -2 . Moreover, the temperature 
pegged at the 200 keV limit, well in excess of a physically 
plausible fit for a CV, indicating that this model is sim- 
ply not a good representation and perhaps a power law 
would be equally good. The "best" (physically plausible) 
fit in this case gave r = 0.5 ± 0.7 with N H = 3± 3 x f 22 
cm~ 2 ; generally poorly constrained and a poor goodness 
of fit of 76% 14 . For high-inclination LMXBs, it is of- 
ten found that there is additional local absorption due 
to material ve rtically extended at the accretion disc rim 
l)Churchll2001|) . while for IPs the material in the accre- 

14 This XSPEC command simulates a user defined number of 
spectra (we used 5000) based on the model and writes out the 
percentage of these simulations with the fit statistic less than that 
for the data. If the observed spectrum was produced by the model 
then this number should be around 50%. This command only works 
if the sole source of variance in the data is counting statistics. 



tion flow can have the same effect. In both cases, the 
absorption may only partially obscure the X-ray source 
spatially and, for the case of an orbitally averaged spec- 
trum as we have here, only for part of an orbit. Hence, 
we included an additional column scaled by its partial 
covering fraction, whilst fixing the other at the nominal 
(maximum) Galactic value. The inclusion of this extra 
component does improve both fits, increasing the pho- 
ton index, but the plasma temperature still pegs at its 
maximum. At the same time, if we fix fcTbr = 30keV 
(typical for IPs) this only increases goodness from 56% 
to 61%. The covering fraction is in fact close to unity 
(0.9 ± 0.1 and 0.95 ± 0.02 respectively). In conclusion, 
our joint fits to the X-ray spectra (as shown in fig.0 and 
detailed in Table. 0J) can only constrain the absorption 
column, requiring a significant local absorption, which 
could equally be accounted for in a mCV interpretation 
or (high inclination) LMXB interpretation. As to the 
emission model, the data seem best fit with a power law, 
but this is far from conclusive. 

5. NATURE OF SDSS J1023 

The observational constraints supplied by our new op- 
tical and especially X-ray data lend significant support 
to the LMXB hypothesis of TA05. The 0.5-10 keV lu- 
minosity of SDSS J1023 is 10 32 ' 4 ergs _1 assuming a dis- 
tance of 2 kpc, which places it squarely in the range of 
luminosities for neutron star LMXBs in quiescence (see 
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Fig. 2.— SDSS J1023: from top to bottom, light curves in 
white light from NOFS, B-band from the XMM-Newton/ OM and 
in the X-ray (0.2-10 keV). They are plotted against phase according 
to our refined ephemeris. The optical curves show smooth, but 
non-sinusoidal profiles, as can be seen by the deviations from the 
overplotted best fit sinusoids. 




l 2 
channel energy (keV) 

Fig. 3.— XMM-Newton/EPIC spectra for SDSS J1023: black- 
pn; dark grey-MOSl; light grey-MOS2. Top: best fit "CV" model 
consisting of a single temperature optically thin plasma model 
(mekal), absorbed by Nh = 1 X 10 9 cm -2 . Bottom: best fit 
"quiescent LMXB" model comprising a soft neutron star atmo- 
sphere (NSA) and a simple power law, with a higher column of 
N H = 4.3 x 10 20 cm" 2 . 



e.g. iCamnana fc Stellal l2fl04D . However, in this energy 
range the flux is domina ted (97%) by the hard power 
law component. Recently. Uonker et al.l l)2004bla|) exam- 
ined the relation between the X-ray luminosity of neu- 
tron star transients in quiescence and the fractional con- 
tribution of the power law component. They found a 
correlation for Lx > 1 — 2 x 10 33 ergs _1 , and an anti- 
correlation for lower luminosities. With its relatively low 
luminosity SDSS J1023 lies at the top end of the anti- 
correlation having a negligible NSA contribution, similar 
to both the LMXB in Terzan 5 and the millisecond pul- 
sar SAX J1808. 4-3658 in quiescence. In the pulsar, the 
leading model for the power law component is the inter- 
action of the pulsar wind with that fro m the donor star 
({Stella et al.lll994HBurderi et alJl200.1) . but as yet there 
is no clear origin for the case of weakly magnetized neu- 
tron stars. Since SDSS J1023 has not been observed with 
a large area X-ray detector, such as the RXTE, during 
an outburst, we have no constraints on its pulsed X-ray 
emission; it is interesting to speculate that it too may 
host a millisecond pulsar. 

Another notable feature of SDSS J 1023 is its short 
timescale (i.e. ~10's of minutes) X-ray variability; again 
this is somewhat unusual, but not without precedent, for 



example Cen X-4 and EX O 1745-248 ijCampana et alJ 
120041 fWiinands eiT al. 2005). As for the spectral param- 
eters themselves: with fcT e g w 50 eV it is lower than 
for most quiescent LMXBs (i.e. 100-300 eV), but fits in 
with the trend for those systems with very weak NSA 
components; the photon index« 1.3 lies well within the 
normal 1-2 range. 

Lastly, now that we have a measure of the quiescent X- 
ray luminosity, we can consider how much brighter it may 
have become during its active state (outburst), and con- 
sider its detectability during previous X-ray observations. 
A search of the HEASARC archives yields no previous 
pointed observations by other X-ray satellites. Naturally, 
its position was covered by both the ROSAT All-Sky Sur- 
vey (RASS) and the RXTE/ ASM, though in neither has 
it been detected. The null result for RASS (upper limit of 
0.03 counts/s) is consistent with our XMM-Newton flux 
measures: we predict merely 0.015 counts/s in the RASS 
0.5-2.5 keV band. However, the fact that SDSS J1023 
remained undetected by the ASM even during its bright 
phase can provide a useful constraint. The nominal ASM 
point source sensitivity is 8.4 x 10~ 10 ergcm~ 2 s _1 , hence 
SDSS J1023 could have brightened by a factor ~2000 and 
still remained undetected by the ASM. Such a brighten- 
ing is well-within the wide spread seen amongst X-ray 
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Fig. 4. — Optical light curves for SDSS J0932. Black filled circles (and error bars) show data from the NOFS 1-m telescope. These data 
were taken in white light; the ~ V indicates that an approximate zero-pointing onto the Johnson system has been applied. Dark grey open 
circles show the B-band data from XMM-Newton/OM.. Seven eclipses are clearly visible, as well as, non-eclipse modulation structure that 
changes significantly between the epochs. 



transients. 

In addition to the X-ray, we can examine its activity 
history via a number of optical observations. Catalog 
data provide snapshots of its state at their epochs of ob- 
servation. We have four independent data points span- 
ning 1952 to 2000, and in each case the approximate B 
and R magnitudes show the source to be quiescent. Fur- 
thermore, all observations since the late 2000/early 2001 
epoch of the Sloan spectroscopy and the Bond et al. se- 
ries of observations show the source to be quiescent. It 
therefore appears that SDSS J1023 has a low-duty cycle 
of activity, which is seen in many other X-ray transients, 
and moreover, fits in with a low-temperature and lumi- 
nosity from the neutron star atmosphere. 

In summary, we can argue that all observations to date 
find a system with characteristics and behaviour fully 
consistent with a LMXB, which spends most of its time 
in quiescence. 

For the CV alternative, we immediately run into prob- 
lems: what provides its X-ray emission in such a qui- 
escent state? A number of CVs are observed to enter 
very low states, optically up to ~7 mags fainter than 
normal: the VY Scl stars. However, they spend the 
majority of their time in active states. Only one IP 



has been observe d to transition into such a very low- 
state, V1223 Sgr l)Garnavich fc Szkodvl 1X988). in a sur- 
vey of the Harvard plate collection, hence nothing more 
is known about its characteristics in that state. Phys- 
ically, in the absence of mass transfer onto a white 
dwarf, the only plausible source of X-rays appears to be 
the chromospheric emission of the fast-rotating donor: 
as is the presumed case for SDSS J155331. 12+551614. 5 
ijSzkodv et a.l .1 12004). But for SDSS J 1023, this option is 
inconsistent with the hardness of its X-ray spectrum. All 
active binaries when fit with optically thin thermal emis- 
sion models yield temperature < 2 keV (see results from 
Chandra observations of 47 Tuc and f rom the RASS: 
iHeinke et al1l2005l IPenTpsev et al.l[l99^) . not > 15 keV 
as we find. Moreover, the ratio of X-ray to optical flux 
is far greater than for any of the RASS or 47 Tuc active 
binaries. 

6. NATURE OF SDSS J0932 

Although the rough V magnitudes of SDSS J0932 (6 
epochs- SDSS photometry through XMM-Newton obser- 
vations) do vary, they differ by much less than a mag- 
nitude, which is likely due to random sampling of the 
eclipses in some of the longer exposures. With no clear 
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evidence of significant state changes (unlike SDSS Jf023) 
a transiently accreting LMXB (or dwarf nova) model 
seems less favored. In any case the high excitation lines 
of its optical spectra suggest either an active LMXB or 
magnetic CV. Furthermore, the high X-ray column (and 
partial covering model) further restrict the options to ei- 
ther a high-inclination LMXB or most likely an IP. As 
we noted, from the X-ray spectral fits alone it is hard to 
choose between these competing interpretations. Hence, 
we have also compared the X-ray to optical flux ratio of 
SDSS J0932 to those of all known IPs 15 and to a sample 
of 8 high inclination LMXBs with periods less than 10 
hrs. Clearly, the exact details of X-ray spectrum, and 
the Galactic absorption are not taken out in these esti- 
mates. Nevertheless, if we normalize our values such that 
SDSS J0932 has F x /F opt = 1 we find it fits well within 
the scatter of values for IPs, which range from 0.1 to 
10, though clustered in the 1-2 range. In contrast, most 
of the high inclination (actively accreting) LMXBs have 
Fx /F opt ~ 100 - 1000; low for LMXBs because the disk 
obscures much of the direct X-ray flux, but still orders 
of magnitude higher. 

The absence of a firm detection of additional period- 
icities in the optical light curve leaves the IP candidacy 
of SDSS J0932 very much a matter of debate. There is 
also our tentative detection of a photometric variation a 
few percent longer than the orbital period, which would 
be difficult to explain in an IP. Much more extensive op- 
tical photometric monitoring is without doubt required 
to resolve these two issues. 

15 See Koji Mukai's compendium, http://lhcawww.gsfc 
. nasa. gov / users / mukai / iphome / iphome.html 



7. CONCLUSIONS 

We have obtained XMM-Newton X-ray and additional 
optical followup observations of two accreting bina- 
ries recently discovered in the SDSS: SDSS J1023 and 
SDSS J0932. Consideration of their optical spectra and 
magnitudes (and history thereof) suggests that we ob- 
served SDSS J1023 during its (normal) inactive state, 
while SDSS J0932 was actively accreting at the time. 

For SDSS J1023 th e X-ray results favor the (transient ) 
LMXB hypothesis of iThorstensen fc Armstrond |2005), 
both in terms of X-ray luminosity, spectral parameters 
and the relative contributions of the best fit neutron 
star atmosphere + power law model. The distance in- 
dependent X-ray-to-optical flux ratio and the estimated 
range of X-ray luminosities are inconsistent with what 
we would expect from an inactive CV, i.e. the chromo- 
spheric emission from the low-mass donor star. 

Our two epochs of time-series optical photometry con- 
firm that SDSS J0932 is indeed a high inclination system, 
exhibiting deep eclipses on a period of 1.57 hr. Unfortu- 
nately, given the very low X-ray count rate of SDSS J0932 
we are unable to construct a useful X-ray light curve. 
Neither are we able to constrain the emission model well 
but our X-ray spectra clearly require significant local ab- 
sorption. This could be accounted for in cither an Inter- 
mediate Polar or high inclination LMXB model. We find 
that its X-ray-to-optical flux is fully consistent with the 
known IPs, but at least an order of magnitude less than 
that of the X-ray faintest LMXB. 
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